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The first aim of this paper is to derive under minimal hypotheses Ehrenfest's theorem concerning the time 
evolution of the Wick symbol for a quantum observable (i.e, its average taken on coherent states). This 
theorem (see [CTDL]) says that the Wick symbol follows the classical mechanics equations as the semiclassical 
parameter h goes to 0. This result shall be precisely stated in section 1 and proved in section 6 for observables 
that are not pseudo-differential operators and consequently in the case where this theorem cannot be deduced 
from Egorov's theorem. One may see [R2], [CF], [CRl], [CR2], [S], [Ul], [U2] for references concerning 
coherent states and some of their applications and Ammari-Nier [AN1][AN2] in the case of the infinite 
dimension. 

Following our preceding work [AKN] , the second purpose here is to study the time evolution of the Wick 
symbol of a trace class operator solution to the time dependent Hartree Fock equation (TDHF). The work 
in [AKN] is devoted to the Weyl symbol of such of an operator ph (t) under the hypothesis that the initial 
operator ph{0) is a pseudo-differential trace class operator belonging to the the class of operators studied by 
C. Rondeaux [R] (actually, the class of operators in [R] is independent on h and the parameter h is inserted 
here in order to consider semi-classical operators). Under weaker hypotheses than those assumed in [AKN] 
we shall study here the Wick symbol of such a solution, instead of its Weyl symbol. A solution to (TDHF) is 
a trace class operator and its Wick symbol is a function belonging to L^(IR^"). We shall prove that the time 
evolution of the Wick symbol for a solution to the (TDHF) equation "tends to the time evolution associated 
to the Vlasov equation" as the limit h goes to 0. This limit shall be understood in the L^(1R^") sense. The 
precise statement of this result is written in section 2 (see theorem 2.1) and its proof is given in section 7. 

One of the tools used in order to established both of these two results is the approximation of a bounded 
operator in L^(IR"), (resp. of a trace class operator), verifying rather weak hypotheses, by pseudo-differential 
operators lying in the class of Calderon-Vaillancourt [CV] (resp. in the class of C. Rondeaux [R]). This process 
is similar to a convolution where translations are now replaced by an action in the Heisenberg group. The 
details are found in section 5. 



In order to motivate the introduction of the Wick symbol it is given in section 3 an explicit example of a 
trace class operator having a Weyl symbol which is not in L^(IR^"). On the contrary, the Wick symbol of 
a trace class operator is always in L^(1R^") (this point is proved by C. Rondeaux [R]) together with all of 
its derivatives (this fact is derived in section 4). Moreover, an asymptotic expansion of the Wick symbol 
for the product of two operators when one of them is a semi-classical pseudo-differential operator and when 
the other one is any trace class operator is written in section 8. The bound for the remaining term of this 
asymptotic expansion is effectuated in the L^(IR^") norm. Note that an analysis for the commutator of an 
operator defined by the Weyl calculus with another operator defined with the anti-Wick formalism is found 
in R. Schubert [S] (theorem 4.1.12). The results in section 8 below can be formally seen as the dual in some 
sense of those in R. Schubert. One may also see lemma 2.4.6 of N.Lerner [LI]. In section 9 it is obtained 
an equation being approximatively satisfied by the Wick symbol of a solution to (TDHF) with an arbitrary 
small error term (see [DLERS] for the idea of such an equation). Finally in section 10, it is observed that the 
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convergence's result established in corollary 2.2 is not uniform on IR even if the two potentials are entirely 
vanishing. 

1. Statement of the result : the case of bounded operators and Ehrenfest's theorem. 

In this section the time evolution of the system is associated to the following Hamiltonian (which depend on 
h>0) 



(1.1) Hh = -h^A + V 

where F is a C°° real- valued function on IR" which is bounded together with all of its derivatives. We 
also denote by Hh the unique self-adjoint extension of this operator. Let {Afi)h>o be a family of bounded 
self-adjoint operators in L^(IR"). The operator Ah{t) corresponding to the evolution of the operator A at 
the time t is 

(1.2) Ah{t)=e'i"^Ahe-'i"'^. 

According to the standard statements of Ehrenfest's theorem (see [CTDL]) the average of Ah{t) taken on 
coherent states is supposed to approximatively follow the classical mechanics's laws as h tends to . 

The coherent state centered at the point X = (x,^) in K^" is the following function, depending on /i > 

(1.3) *x,ft(«) = (7r/i)-"/^e-^Tet«-«-*^-« X = (x,^) e IR^". 
The average taken on ^xh of a bounded operator AmH = i^(IR") is 

(1.4) ar'^iAKX) =< A^xh, ^xh > ■ 

The function af^'^'^^{A) is called the Wick symbol of A. When A is a bounded operator we shall notice 
(see proposition 4.1) that the function a^^'^'^^A) is C°° on IR^" and we shall give precise estimations on its 
derivatives in terms of the parameter h and of the norm of A. Let us now specify in which sense does this 
function follows the classical mechanics equations under the limit h tends to 0. 

Let i^t{x-, C) = {(lt{x-i C)^Pt{x, ^)) denotes the Hamiltonian flow of the function H{x, ^) = |^|^ -|- V{x) starting 
at [x,^), i.e., the solution to, 

(1.5) 4{x, = Mt, X, Ptix, = -^Viqit, X, 0) 

such that qo{x,£,) = x, po{x,^) = ^. One can say that a function f{x,^,t) "follows the classical mechanics 
laws" when f{X, t) = fi^X), 0) for all X e IR^" and all t e K. 

If the initial data Ah is not assumed to be a semiclassical pseudo-differential operator (see D. Robert [Rl]), 
then the function a'^^'^'^ {Ah{t)) may not have a limit as /i ^ 0. Nevertheless, we shall precisely state in which 
sense this function "follows the classical mechanics's equations under the limit h tends to 0". Our goal is to 
prove under rather weak hypotheses and for all time f S IR that, 

^wickf A „wick / 



(1.6) lim ar^^AhmX) - ar%AH{0)){MX)) 



= 

^0 ■ ■ " 

In particular, this is not fT^*'^'^(A/j(f))(X) that has a limit when /i — ^ but the difference between this function 
and another one which follows the classical mechanics laws and takes the same value as a'^^'^'' {Ah,{t)){X) at 
time 0. 
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Let us now precisely state the hypotheses implying the validity of the limit (1.6). We shall denote by 
W"J'(]R^") the space of functions in LP(]R^") having all of its derivatives up to order m in LP{M'^") (1 < 
p < +00, < m < +00). 

First, let us recall how the Egorov theorem is classically applied the case of operators Ah which are semiclas- 
sical pseudodifferential operators. In this case, Ah is associated, through the Weyl calculus, to a function 
Fh in W°°°°{m?''), which is bounded in T^°°°°(]R2n) independently of h G (0, 1]. Namely, for all feH, 

(1-7) {Ahf){x) = {2nh)- [ Fh{^^,Oe^^^-y^-^f{y)dyd^ 

According to C alder on- Vaillancourt [CV], Ah is bounded in H. We agree to express (1.7) under the two 
following expressions 

(1.8) Ah = Opl'y\Fh) Fh = a^'^Ah) 

(The Weyl symbol of an arbitrary bounded operator being a priori in H is a tempered distribution on M^".) 
If Fh is either independent of h or admitting an asymptotic expansion in the powers of h as h tends to 0, 
Egorov's theorem shows that: 

(1.9) lim (ar^AhimX) - cr^' (AhmiMX))) = 
In that case, the limit (1.6) then follows. Indeed, 

(1.10) ar''{Ah{t))=ei^ar\Ah{t)) 

(Even if this result is standard, we shall nevertheless give a proof in proposition 4.4.) Besides, Egorov's 
theorem also shows that, if Ah is written as in (1.7) with Fh in W°°°°{JR^'^) and uniformly bounded in 
h e (0,1] then a^^^^^Ait)) is also a function in W°°°°{JR^''), bounded uniformly on h G (0,1]. As a 
consequence, 

(1.11) hm (ar^AhitrnX) - (7r^'(A,(i))(X)) = 
for all t e ]R. In this case, the limit (1.6) then follows from (1.9) and (1.11). 

The assumption A = Op'^''^'{F) with F in VF°^°^(]R2») is rather strong and may be expressed in terms of 
commutators in a standard way. Let Pj (h) and Qj (h) be the momemtum and the position operators, 

(1.12) p.(/,) = ^_^ Qj{h) = xj 

According to Beals characterization result [B] , the fact that Ah may be written as Ah = Op^^^^ (F) with F 
in W^°°°°(]R^") is equivalent to the fact that the commutators {ad P{h))°'{ad Q{h)YA are bounded (for all 
multi-indexes oi,j3). Saying that F (which may depend on h) stays bounded in W°°°°{^") is equivalent to 
the fact, 

h-^\^\+\^^{ad P{h)r{ad Q{h)fA\\c(H) < 

with Majj independent on h. 

We are first concerned with the proof of the limit (1.6), which may possibly be viewed as a form of Ehrenfest's 
theorem, but under much weaker hypotheses than in the case where Ah is a semiclassical operator. Namely, 
only single commutators (instead of iterated commutators) of the operator Ah with the operators Pj [h) and 
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Qj{h) defined in (1.12) are assumed to be bounded operators. Our estimates shall involve the following 
expression, 

n 

(1-13) iriA) = + 



The theorem below provides the inequality on which rely the proof of Ehrenfest's theorem with weakened 
hypotheses. 

Theorem 1.1 For all operators A in such that the commutators [Pj{h)^A\ and [Qj{h),A\ are hounded 

operators inT-L (1 < j < n), the operator Ah{t) defined in (1.2) and the Hamiltonian flow tpt defined in (1.5) 
satisfy 

(1.14) \ar\An{t)) - {ar\A)) o ^t\\^^^^^^^ < C{t)^I^{A) 

where 1 1-> C{t) is a function defined on IR, hounded over any compact set, depending only on n and on V 
and where I^^iA) is defined in (1.13). 

In order that (1.6) can be viewed as Ehrenfest's theorem it suffices to replace A by a family of operators 
Ah such that the right hand-side of (1.13) tends to as h goes to 0. This assumption is satisfied when 
Ah = Opl^'"\Fh) with Fh in W°°~(]R2") and uniformly bounded in /i e (0, 1]. 



2. Statement of the result: the case of trace class operators and the time dependent Hartree 
Fock equation. 

Before introducing the time dependent Hartree Fock equation (TDHF) let us first speciiy the space to which 
the solution belongs. We shall denote by C^(H) the space of all trace class operators onH = L^(]R"). We 
shall denote by V the subspace of operators A in C^CH) such that the commutator [A, A] is also in jC^{H) 
where A is the Laplacian. 

We consider two real- valued functions V and W in M^°°°°(1R"). For all /i > we shall say that a function 
t — > ph{t) being from IR into C^{H) is a classical solution to the time dependent Hartree Fock (according 
to the terminology of Bove da Prato Fano [BdFl,BdF2]) when this mapping is also continuous on IR into X> 
and if, 

(2.1) ^h^^Ph{t) = -h^[A,ph{t)] + [V^{ph{t)),Ph{t)] 
where Vq{ph{t)) denotes the multiplication operator by the function: 

(2.2) Vq{x,ph{t)) = V{x)+Tr{W^ph{t)) 
where Wx is the multiplication operator by the function y W{x — y). 

We consider a family {ph{t))h>o of classical solutions to (TDHF). We suppose that the operator ph{0) is 
trace class, self-adjoint > with a trace equal to 1. We set: 

(2.3) UhiX, t) = {2nh)-^arHphmX) 
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As we shall see in section 4 this function lies in W°° {JR it is > and 

(2.4) / Uh{X,t)dX = l 

We denote by Vh{X,t) the solution to the Vlasov equation, 

/r, r.\ ^^'^ ^ o t ^^f" \^ dVcl{x,Vh{.,t)) dVh _ 

such that 

(2.6) Vh{X,0) = Uh{X,0) 
We have set, 

(2.7) Vci{x,Vh{.,t)) = V{x)+ [ W{x-y)vhiy,v,t)dydr] 
The function Vh{-,t) is itself in L^(]R^"), it is > and its integral equals to 1. 

The counterpart of Ehrcnfest's theorem for a family of solutions to (TDHF) consists into saying that the 
function Uh{-,t) defined in (2.3) satisfies Vlasov equation "under the limit h tends to 0". This point of view 
may specified in diff'erent ways. In this section our aim is to compare as in theorem 1.1 the functions Uh{-,t) 
and Vh{.,t) and to show that, under suitable assumptions on p?i(0), their difference in norm tends to 0. Since 
we are concerned with trace class operators it makes sense that the norm involved in the estimation of the 
diff'erence Uh{-,t) — Vh{-,t) is the L^(]R^") norm. A first answer to this issue is given in our preceding article 
(preprint) [AKN] where we assume that the operator ph(0) is a pseudo-differential operator belonging to 
a class of operators introduced by C. Rondeaux [R] (the only modification here being the insertion of the 
parameter h). The article [AKN] also gives an asymptotic expansion at any order of Uh{-,t) — Vh{-,t). We 
consider here this problem with a weaker hypothesis than the one in [AKN] . We arc only assuming that all 
the commutators [Pj{h), ph{0)] and [Qj{h), Ph{(i)] are trace class operators. All these estimates shall use the 
expression, 

1 " 

(2.8) itiPhiO)) = -f^J2\\^Pj{h),Phm\\cHn) + \\[Qj{h), PhiOMcHn) 

Theorem 2.1. Let {ph{t))h>o be a family of classical solutions to (TDHF) with real-valued potentials V and 
W in H^°°°°(]R^"). We suppose that the operator Ph{0) is trace class, self-adjoint > 0, with a trace equal 

to 1. We assume that all the commutators [Pj{h), ph{0)] and [Qj{h),pfi{0)] are trace class operators. Then, 
there exists a function t i— >■ C(t), bounded on any compact set o/IR such that, for all h G (0, 1] and for all 
t e JR, 

(2.9) \\uh{.,t)-Vh{.,t)\\L^(m^) < C(t)V^/r(p^.(0))e^W^^'('"'(°» 

Corollary 2.2. Under the assumptions of theorem 2.1, if Ifl'{ph{Q)) remains bounded when h tends to 0, 
then we have, 

(2.10) lim ||Mft(.,t)-Uh(.,t)||ii(]R2„) =0 
for all t e H. 
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If phiO) is a pseudo-differential operator written as ph{0) = {2Trh)"Op,"''" (Fh) with Fh in I4^°°1(IR2") and 
uniformly bounded in h then the assumption in corollary 2.2 is always satisfied according to the analogous 
result of the Beals's characterization given in [R] (see [AKN] for dependence on the parameter h). 

Theorem 2.1 shall be proved in section 7. An other way to give an insight into Ehrenfest's theorem consists 
into proving that the function Uh{X,t) satisfies an equation that formally shrink to the Vlasov equation 
when h tends to 0. It is this point of view that we shall develop in section 9. 



3. A counter-example. 



By making explicit an idea of C. Rondeaux [R] we shall give an example of a trace class operator having a 
Weyl symbol not being in L^(1R^"). Let us mention before some properties for the Weyl symbol of a trace 
class operator that lead us to the choice of the counter-example. 

Proposition 3.1. If A is trace class, then its Weyl symbol is a continuous function on IR^" going to at 
infinity and also belonging in L^(]R^"). If this function is also > then it is necessarily in L^(]R^"). 

Proof. If A is trace class then it is also Hilbert-Schmidt and it is well-known that its Weyl symbol belongs 
to L^(]R^") (cf D. Robert [Rl]). This symbol is bounded since when A is trace class we have, 

(3.1) ar^'(A)(X) = 2"Tr(^ o l^xh) X e IR^" 

The fact that it is continuous and that it is going to at infinity is easily verified using (3.1) for an operator 
written as m — > A{u) =< u,ip > ip, where ip and ip arc in (S(IR"). One then concludes this point by density, 
first for any finite rank A, then for any trace class operator since the set of finite rank operators is dense in 
C^{n). For all functions F in H^°°1(1R2") and G in H^°°°°(IR2") it is well-known (cf D. Robert [Rl]) that, 

(3.2) TriOp'l^'y^iF) o Op™*^'(G)) = (27r/i)-" / F{X) G{X) dX 

Let A be a trace class operator and G be a C°° function on IR^" with compact support. From theorem 
5.2 there exists a sequence of functions (Fj) in W°°^(IR^") such that the sequence of operators Op^^^^Fj) 
converges to A in C^{H), implying from (3.1) that Fj converges uniformly to crJ^*^'(A). It is then deduced 

that, 

(3.3) Tr{A o Opl^y\G)) = (27r/i)-" / al''y\A){X) G{X) dX 

We replace G by an increasing sequence Gjv of G°° functions on IR^" with compact support and converging 
pointwise to 1 when A'' goes to -|-oo, the functions |9"9^Gjv| being all uniformly bounded on N. When 

(tZ''^\A) > we deduce from (3.3) that, 

< {2nh)- j ar^'(A)(A) G^(A) dX < \\A\U^^n) \\Opr'iG^)\\cm 

According to Calderon-Vaillancourt the right hand-side remains bounded as N tends -|-oo. If the function 
a^'^y^iA) is > then it is therefore in L1(IR2"). 

□ 



Proposition 3.1 and the results of C. Rondeaux [R] suggest the construction of an example of a trace class 
operator having a Weyl symbol not being in L^(IR^"). We may assume that h= 1. Let a be a real number 
such that ^ < a < n. We define the operator P by: 



(3.4) P = Opr'(p) p(^,^) = __^___ (x,OeIR^" 
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The function p is not in L^(IR^"). For all A > we define an operator A\ by 
We have: 

1 /"^ 

(3.5) P = -— / e-^X'^-'AxdX 

r(a) Jo 

provided that we verify the convergence. From definition (1.7) the integral kernel of A\ is the function K\ 
defined by: 



Kx{x,y) = (27^)-"^^ e'(--^)-«aA (^^c) d^- 



An explicit computation shows that: 

Kx{x, y) = (27rA)-"/2 e-^a(X)\A'+b(X)\y\-+2c{X)..y) = 7 + T KA) = c(A) = ^. 

4 A 4 

Wc can express the operator A\ as a product, 

(3.6) Ax = (27rA6(A))-"/2 r„(;,) o o T-^^ o S 

where for all a > 0, {Taf){x) = f{xs/a), Sf{x) = f{-x) and: 



-L 



c(A) 



(3.7) {Bxf){x) = I e-(l-l'+l^l')+2''W-V(y)rf2/ m(A) - 

iR" Va(A)6(A) 

The operator B\ is self-adjoint > (since < /i(A) < 1 and B{X) may be identified, up to a multiplicative 
constant, to the exponential of an harmonic oscillator). We then see: 

From (3.6) we have, 

WAxWc^^n) < (27rA6(A))-"/2 ||r„(;,)||£(^) \\Bx\\c^^n) IK^lM^y 

The right hand-side is polynomially increasing as A +00 and it is ©(A""/^) when A tends to 0. With the 
choice of a > -j, the integral (3.5) converges in norm in C^(W) and it properly defines a trace class operator 
P, the Weyl symbol of which is the function p defined in (3.4), and is not belonging to L^(IR^"). 

For trace class operators, it is natural to give a variant of Ehrenfest's theorem where the limit is understood in 
the Z/^(1R^") sense. The above counter-example shows that the Weyl symbol is not suitable for this purpose 
without any supplementary hypothesis. However, proposition 4.2 below shows that the Wick symbol of a 
trace class operator is C°° from IR-^" into L^(IR-^") together with all of its derivatives. 

4. Differentiability of the Wick symbol. 

Proposition 4.1. Let A he a hounded operator in H. Then the Wick symhol of A, namely, the function 
^wick^j^^ dgjljigfl ]p>^2n Cj^^j^ ig a C°° function on IR^", hounded together with all of its derivatives. 
For each multi-index a there exists Ca such that: 

(4.1) II Wr'(^)l|L~(ni-) < a/i-l"'/^ WAWany 
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Proof. We shall use the following function, 

(4.2) 5.(A)(x,y)= <f^'-/^'-> , 

defined on IR^" x M^", where the functions ^xh are given in (1.3). Direct computations shows that: 

(4.3) < ^xh, ^Yh >= e-ikl^-^P+AM^.7)_ 
In particular, 

(4.4) I < ^xh,^Yh > I = e-TCl^-'>"l' ll^-^^ll = 1. 
Consequently: 

(4.5) \Sh{A){XX)\ = e^l^-^l'l < A^xh,^Yh > \ < e^^'^-^l' \\A\\ciny 
An important property verified by coherent states is that: 

(4.6) <f,9>= (27r/i)-" / < /, >< 5 > dX 
for all / and g inH. Applying this equality, we see that: 

(4.7) ar'\A){X) = Sh{A){X, X) = {2Trh)-'" [ Bh{X, U, V, X)Sh{A){U, V)dUdV 
where for all {X, U, V, Y) we have set: 

(4.8) B,iX, U, V, Y) = <^xH,^YH> • 

Prom (4.3) we have: 

(4.9) Bh{X,U,V,X) = e-2t(^-n(^-£7)-^|t/-v'l^ 
We verify that: 

where 

K^{X, U,V,h) = (^l + \^-Ill±^i^^ ^ e-i-^i\x-u\^Mv-x\^). 

Consequently: 

(4.10) V^Vr"'=(A)(X) < C„(27r/i)-2" / I < A<iuh, "^vh > \Km{X, U, V, h) dUdV. 
By bounding from above | < A^uh, ^vh > \ by the norm of A we then obtain (4.1). 

Proposition 4.2. If A is in C^CH), then its Wick symbol ah{A) belongs to W°°^ {JR?'^) . For each multi-index 
a there exists Ca > such that, 

(4.11) (27r/i)-"|| V™ar'=(^)IUi(iR-) < Cm/^"^ WMcHny 
Inequality (4.11) is proved when m = by C. Rondeaux [R]. 



Proof. Similarly to proposition 4.1, we have inequality (4.10). Integrating with respect to X we obtain: 



G™([/) = (1 + 1 (71)^6- 

Since the function Gm is in L^(]R^"), proposition 3.2 is a consequence of the following lemma. 
Lemma 4.3. iet A be a trace class operator and G be a function in L^(IR^"). Then we have: 

(4.12) (27r/i)-2»^^^ I < A^xh,'i!Yh > G (^^) | dXdY < (27r)-"||G|Ui(R2n)||A||£i(H). 



Proof. We can write A = B1B2 where i?i and B2 are two Hilbert- Schmidt operators. Using the fundamental 
indentity (4.6) verified by the coherent states we see that, for all X and Y in K^", 

< A^xh, ^Yh >=< B2^xh, B^^Yh >= (27r/i)-" / uzh{X) vzh{Y) dZ 

where we have set uzh{X) =< 32^^ xh, zh > and vzh{X) =< "^zh, -B^^I/xft >• Let 1^ be the left hand-side 
of (4.12). From Schur's lemma: 

h < (27r/i)-3"/i-||G|Ui(i,2„) / WuzhWmiR^r^) \\vzh\\mM--)dZ. 

From (4.6), we have \\uzh\\L^{TR-^) = (27r/i)"/2||S$*2^|| and Wvzhh^M-^) = (27r/i)"/2||Si*zft||. Conse- 
quently, 

h < (27r/i)-2"/i"||G|Ui(K2n) / \\B*2^zh\\ dZ. 

From the basic indentity (4.6) for coherent states, 

(27r/i)-" / WBj^HzhfdZ = {2nh)-'' [ < B* Bj^H zh,^ zh > dZ 

= Ti{BjBj) = ||-Bj|||2(^) 
where ||Bj||£2(9^) is the Hilbert-Schmidt norm of Bj. Therefore, 

Ih < (27r)~"||G||ii(iR2„) ||Bi||£2(^) ||B2||£2(^). 

Taking the infimum over all A written as ^ = -B1-B2 we then obtain (4.12). 



□ 



Let us finally recall the following proposition even if it is standard. 

Proposition 4.4. The Wick symbol a^^^''{A) of an operator A is related to its Weyl symbol a^^^\A) by 
(4.13) ar'\A) = ei^a;;'y\A) 

where A is the Laplacian on H^". 
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Indeed, setting Fh = a^^^\A), the expression in definition (1.7) for the Weyl calculus may be written as: 

(4.14) A = (tt/i)-" / Fh{Y)^YhdY, 

where, for all Y = {y,r]) in IR^", T,Yh is the operator ("symmetry operator") defined by, 

(4.15) {^Yhf)iu) = e^(«-y)«/(2y - u) Y = {y, r,) G K^". 
A direct computation shows that: 

(4.16) af^'=(Sy^)(X) =< Syfcfx/., ^xh >= e-^^. 
Equality (4.13) then follows. 

5. Pseudo-difFerential approximation. 

The class of operators having a Weyl symbol in W^°°^(]R^") (introduced by C. Rondcaux) is dense in the 
set of trace class operators similarly to the fact that M^°°^(]R") is dense in L^(IR"). The approximation 
process has a strong connection with convolutions and may probably be also applied to Schattcn's classes. 
This process is often employed in representation theory. However, as for a function in L°°(]R") needs to be 
continuous in order to be in the closure of W°°°° {WC^) , a bounded operator needs additional hypotheses in 
order to be in the closure of the class of Calderon-Vaillancourt for operators. 

For all X = {x, ^) in ]R" and for all /i > let Wx,(,h be the operator defined by 

(5.1) {Wx,hf) = {W,,i,hf){u) = f{u - x)ei^-^-^^-i 

for all / G i^(]R"). It is a common representation of the group of Heisenberg. Thus, the coherent state ^xh 

verifies, 

(5.2) ^x,h = Wx,h'^o,h 
We have, for all X and Y in IR^n 

(5.3) Wx,hWY,h = e^^^'^'^^Wx+Y^h 

where a is the symplectic form a{{x, ^), {y, t])) = y.^ — x.rj. For all operators A in £('H) and for all /i > let 
us define, 

(5.4) TftA = (tt/i)-" / e-'-^Wx,hAW*x,hdX 

We begin with the case of bounded and trace class operators satisfying some hypotheses involving their 
commutators with position and momentum operators. 

Theorem 5.1. a) We have: 

(5.5) WThAWau) < man) 
for all operators A in C{'H) and for all h > 0. 

b) When the commutators [Pj{h),A\ and [Qj(h),A\ are hounded operators we have: 

„ n 

(5.6) ||A-r,A||£(„)<-^^ll[p,(/i),A]|U(^) + ll[0,(/i),A]||^(^). 
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c) When the operator A is trace class together with the commutators [Pj{h),A] and [Qj{h),A], we have: 

„ n 

(5.7) WA-ThAWc^u) < + 

d) The Wick symbol of the operators A and ThA are related with: 

(5.8) ar'=(rfc^) = e^^ar'=(A). 
The Weyl symbol of ThA is equal to the Wick symbol of A. 

Proof. Point a) is clear since Wxh is unitary. For any 9 in [0, 1] define: 

T{e, h)A = (tt/i)-" / e--^W0x,hAWex h dX. 
Thus T(l, h)A = ThA and T(0, h)A = A. We verify that: 
d 

Wex,hAW0x,h = ^ 

Consequently: 



g^Wex,hAW^x^h = lT.[^iWgxAPM,A]We\h- ^jWexAQj{h),A]We\h 



P-r,A||^(„)<if^(^/i)-" / e--^\\x,\ mih),A]\\cin) + \^A \\mh),A]\\cin) 



dxd^dO 



C 



^ 7r E \\{(^dPj{h))A\U^n) + \\{adQjih))A\\c^n) 
proving point b) and also point c) with straightforward modifications. For the point d) we see that: 
ar^''{rhA){X) =< (r^^)*xft,*xft >=< {ThA)Wx,h'ifoh,Wx,h'ifoh > 

= (tt/i)-" / < WY,hAW^hWx,h'i^oh,Wx,h^oh > dY 

= (tt/i)-" / e"^ < AWx-Y,h'i^oh,Wx-Y,h'i^oh > dY. 
We have used here (5.3). Consequently, 

a^'^''{ThA){X) = (tt/i)-" / e-^<T^''=^(A)(X - Y)dY 

which is (5.8). According to proposition 4.4 we also have: 

ar'^{ThA) = ei^aj:'^\ThA). 
Since the operator e^'^ is one to one we then deduce as it is mentioned, 

weyl/^ A\ wick / a\ 

(^h i'hA) = (Jf^ [A). 
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Next wo consider the case of trace class operators without additional assumptions. The result below does 
not have any counterpart in the case of bounded operators. 

Theorem 5.2. The space of operators written as 0P^^^\F) with F in W°°^{JR^"') is dense in the space 
C^{T-L) of trace class operators. 

Proof. For this purpose, we modify the approximation process and we set: 



e — —Wx,iAW*xi dX 



for all A > and for all trace class operators A. Let us show that: 
(5.9) lim||r,'(A)-A||^i(«)=0 

for all trace class operators A. Since we clearly have 



A = {-kX)- 



e — —AdX 



IR2" 



then we see that: 



\\Ti{A) - < (ttA)-" / e-— \\Wx,hAW*x,h " ^Ikn^) 

J\X\<5 



... + ( 



ttA)-" / 
J\x 



e ^ 



\x\>s 



dX 



\Wx,xAW;c,h\\cHn) + \\A\\cHn) 
for all 5 > and for all A > 0. For all trace class operators A and for all £ > there exists S > such that: 

\X\<S,^ \\Wx,iAWl,-A\\ciin)<e. 

Indeed, this property is first verified when A is of the form f ^< f,ip > tjj with ip and ip in <S(]R"), it is next 
derived by density for finite-rank operators and then, by density again for trace class operators. Besides, 
d > being fixed, we have: 



lim ( 



ttA)-" / 
J\x 



\x\>s 



\Wx,iAW^^^\\ciin) + \\A\\cHn) 



dX = 0. 



The limit in (5.9) is then easily obtained. From (5.3), for all X in K^" we have: 



Wx,iTl{A)W^ 1 = (ttA)-" / e-^Wx+Y,iAW^+Yi dY 



(ttA)- 



e--^^Wz,iAW*zi dZ. 



Consequently, for all A > 0, the mapping X VK^.i^K^)^^'!,! is C"^ from H^" into C\H). We see, setting 
X = (x,^) as the variable of R^", 

A.Wx,in{A)W^^, = -Wx,i[Pj{l),n{A)]w*^^, 
A.Wx,in{A)Wl^ = Wx,i [Qjil),niA)\wl, 
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where -Pj(l) is the operator of differentiation with respect to Uj and Qj{l) is the multipheation operator by 
Uj. Consequently, all order iterated commutators Tl{A) with the position and momentum operators Pj{l) 
and Qj(l) are trace class. Prom the result of characterization of C. Rondeaux [R] (the analogous one of 
Beals's characterization for trace class operators and recalled in the first part of this work) it follows that 
Tl{A) is written as Opr^'(^A) with Fx in W^^im?"-) for all A > 0. 

6. Proof of theorem 1.1. 

Let A be an operator in JC-CH) satisfying the assumptions in theorem 1.1. Set Fh{X) = (t^*°'^(^)(X). 
This function is C°° from proposition 4.1. Denoting by cpt the Hamiltonian flow associated to the function 
H{x,^) = + V{x), we shall use the following function and operators: 

(6.1) WhiX,t) = FhiMX)) Chit) = OpZy\wh{.,t)). 
We shall also use the operator ThA appearing in section 5 and the following operator: 

(6.2) Bh{t) = e'i"-rh{A)e-'i"^. 

We shall compare Wick symbols of the operators Ahit), Bh{t) and Ch{t) and then compare the Wick symbol 
of C{t, h) with the function Wh{-,t). This is accomplished in the three steps below. 

First step. We have, from theorem 5.1, 

(6.3) hr^Ahit) - S,(t))|U~(K-) < WAuit) - Bh{t)\\c(,n) = \\A - T,(^)||£(„) 

... < C^rhI^{A) 

where I^{A) is defined in (1.13). 

Second step. The comparison of Bh{t) and Ch{t) comes from Egorov's theorem. Nevertheless this requires 
some precisions due to unusual estimates satisfied by the derivatives oi Wh{.,t) that we first need to specify. 
Since w/i(.,0) = cr^''''=(A) we deduce that: 

Applying proposition 4.1 to the above commutators we see when fc > 1 that: 

iiv'=u;,(.,o)iUo.(H2.) < a^/ii-('=+i)/2/r(A). 

The derivatives of order > 1 of the Hamiltonian flow cpt associated to the symbol H{x,£^) = + V{x) are 
bounded in i?^" with a bound equaling to 0(1 + i^). Then there exists Mk{t) such that: 

(6.4) ||V'=«;h(.,t)|U=o(n,2„) < Mk{t)h^-^''+^^l^ I^{A). 
The operators Bh{t) and Ch{t) satisfy: 

(6.5) -ih^Ml = [Hh,B{t,h)\ 

-ih^^ = -ihOpZy\dtWH{.,t)) = -ihOpZy\{H,wn{.m- 

From a standard result on the Weyl calculus recalled in proposition 3.2 of [AKN] (first part of this work), 
for all functions F and G in W°°°°(]R2n), the operator R^j^\f, G) defined by: 

[Opr\F), Opr\F)] = ^Opr\{F, GD + R<j^\F, G) 
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satisfies the following estimate: 

(6.6) 0^\F,G)\Uin) <C ^ V^F|Uo.(h..) \\y'Gh^^j,,. 



j>2,k>2 
4<j+k<6n+S 



With these notations, one may write: 

(6.7) -'^^^ - l^h,Ch{t)] = R^^\H,WH{.,t)). 
Note that i?)^ ' {A,Wh{-,t)) = and consequently 

We then may apply inequality (6.6) with the functions F = V and G = Whi-jt). The inequality (6.6), those 
in (6.4) which are verified by Wh{-,t), and the fact that all derivatives of V are bounded allows us to write: 

(6.8) \\R'^\H,Gh{.,mcin) < M{t) h"^ I^{A). 

Prom theorem 5.1, the operator ThA appearing in section 5 has a Weyl symbol equal to the Wick symbol 
of A. Consequently, the Weyl symbol of BhiO) = Th{A) and the one of C?j(0) which is Fh = ct^*'='=(^) are 
equal. Thus, 

(6.9) 5^(0) = C/,(0). 
Prom (6.5), (6.7) and (6.9), Duhamel's principle implies: 

ih[Bh{t) - Ch{t)] = Te^^^" R'j^\H,Gh{.,t)) e"*^^" ds. 
Jo 

Consequently, when t > 0: 

\\Bh{t) - Chimcin) < lJ^0^\H,Gh{.,s))\\cin)ds. 

We then deduce that: 

(6.10) War^^B^it) - Cft(i))lUoo(R2») < WBuit) - Ghimcin) < M{t)Vh I^{A). 
Third step. From proposition 4.4, 

<jr\{CH{t)) = e*^ar^'(C,(f)) = e*^^z;,(.,t). 

Then 

In view of the estimates (6.4) satisfied by Wh{-,t), we obtain: 

(6.11) \\ar\Ch{t)) - «;,.(., i)|Uoo(n,.„) < M{t)VhI^{A). 

Since Wh{-,t) = (^(t^^"'^ {A)j o (p^, then inequality (1.14) in theorem 1.1 arises from (6.3), (6.10) and (6.11). 
7. Proof of theorem 2.1. 
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Let ph{t) be a family of solutions to the equation (TDHF) (2.1), satisfying the hypotheses of theorem 2.1. 
Let Uh{-,t) be the function defined in (2.3). Let Vh{-,t) be the solution to Vlasov equation (2.5) such 
that Vh{-,0) = Uh{-,0)- We shall use the functions Vq{.,ph{t)) and Vci{-,Vh{-,t)) defined in (2.2) and (2.7) 
respectively. We shall also use the following functions 

(7.1) Hnx,^,t) = \^\^ + Vg{x,ph{t)) Hr{x,i,t) = \e + Vci{x,Vh{.,t)) 
and the associated operators through the Weyl calculus, namely: 

(7.2) H^^ = -h^A + V,{.,ph{t)) = -/I'A + V,i{x,VH{.,t)), 

by using the same notation for the function and the corresponding multiplication operator. We shall denote 
by Wh {X, t) the solution to 

(7.3) ^{.,t) = {Hf^^{.,t),WH{.,t)}, 
such that 



(7.4) Wh{.,Q) = Vh{.,Q) = Uh{.,Q). 

In order to compare Vh{-,t) with Wh{-,t) we note that the Vlasov equation (2.5) is written as 

(7.5) ^{.,t) = {Hr{.,t),vU.,t)}. 
We shall use the operator Bfi{t) solution to 

(7.6) ^'^^^ = i^h^i*)' B^{0) = Th{phm, 

where Th is the mapping used in section 5. Finally, we shall also use the following operators: 

(7.7) Ch{t) = {2T,hrOpZy\wh{.,t)) Dh{t) = {27ThrOpZ''ivh{.,t)). 
According to the point d) in theorem 5.1, we have 

(7.8) cTr\B,{0)) = ar^(TMO))) = crrHPhiO)) = {2nhru,{.,0). 
Consequently, 

(7.9) Bh{0) = Ch{0) = Dh{0) 

Theorem 2.1 is a consequence of the comparison between the Wick symbol of the operators Ph{t) and Bh{t), 
between those of Bh{t) and Ch{t), between those of Ch{t) and Dh{t), and finally between the Wick symbol 
of Dh{t) and the function Vh{-, t). Each of these comparisons shall be written using the expression lf^{ph{0)) 
defined in (2.8), and shall correspond to one step of the proof, but before that, we need three more lemmata. 

Lemma 7.1. Let Ph{t) be a family of solutions to the (TDHF) equation satisfying the assumptions in 
theorem 2.1. Let Vh{.,t) be the function defined above. Then, for all integer numbers k> we have 

(7.10) Wy'v^i-MLHM^-) < Ck{t)h-''/'\\pHmcHn) 
and for all integer numbers k> 1 

(7.11) \\v'^vh{.,t)\\LHM^'^) < Ck{t)h-(^-^y^ii^{phm- 
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These estimates remain valid when replacing the function Vh{-,t) by the function Wh{-,t). 



Proof of the lemma. Since we have (7.4) for t = 0, then the estimates (7.10) and (7.11) when t = come from 

proposition 4.2 applied with the operator ph{0) (for (7.10)) and with the commutators [Pj (h) , pfi{0)] and 
[Qj{h),Ph{0)] (for (7.11)). Since the potentials V and W are in VF°°°°(IR"), then the functions Vg{.,ph{t)) 
and Vci{-,Vh{-,t)) are uniformly bounded together with all of their derivatives. Consequently, the estimates 
satisfied at t = by Vh{.,0) = Wh{.,0) remain valid along the Hamiltonian flows associated to the two 
symbols Hj^^ t) and Hj^^{., t). Thus, the estimates (7.10) and (7.11) remain valid for all t for the function 
Vh{-,t) and the function Wh{-,t). 

Lemma 7.2 If W is a function in H^°°°°(IR"), if we denote by Wx the multiplication by y ^ W{x — y) and 
if A is trace class, then 

(7.12) e4^^Tr(W, o A) = (27r/i)-" / W{x - y)ar'''{A){y, r])dydr]. 



Proof of the lemma. From (3.2), if W is as in the lemma and if A = Op^^'^^' (F) with F in 1^~1(IR2"), then 
we have: 

Tr{Wx o A) = (27r/i)-" / Wix - y)aZ''\A){y, r,)dydr, 
ei^-Tr{WxoA) = {2Trhy [ W{x - y)ei'^-a];^^''^{A){y,r])dyd'n. 



and then 



Besides 



/ W{x - y) {ei^n _ /) ar'^\A){y,n)dydv = 0, 



and taking into account proposition 4.4 we deduce (7.12). Suppose now that A is an arbitrarily given 
trace class operator. Then theorem 5.2 shows that there exists a sequence of operators Aj, written as 
Aj = Opl'^y^Fj) with Fj in VK°°1(IR2") converging to A in C^iU), {h > being fixed). Equality (7.12) 
valid for all the Aj is also true for A when taking the limit while using proposition 4.2. □ 



Lemma 7.3. With the above notations, we have: 

(7.13) \\vH{.,t)-wn{.,t)\\L^(.^--) < C{t)lf:{ph{0)) 



h+ |kh(-,s) - Wh(-,S)||L1(IR2")C?S 

J[0,t] 



Proof of the lemma. We deduce from (7.3) and (7.5) that: 

d{vH - «;.) ^ {fff ^(.,t), iv,{.,t) - Wf,i.,t))} + {(Ff ^(.,0 - H):^{.,t)),VHi.,t)} 



dt 

Prom Duhamel's principle, since Vh.{-,Q) — Wh{.,0) = and since the Hamiltonian flow associated to the 
function Hf^^{.,t) preserves the norm of L^(]R^"), we obtain 

\W{.,t)-Wh{.ML-(jR?") < / \\{Hh^{.,s)-Hl'^{.,s) , ?;/,(., s)}||ii(jR2„)ds. 

J[0,i\ 

We have: 

fff ^(., s) - Hj^'^i., s) = {I- ei^=')TriWx o p^{t)) + e^^^TriW, o p^{t)) - [ W{x - y)vh{y, V, t)dydv. 
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From lemma 7.2, we have: 

e^'^-TriW^ o ph{t)) = I W{x - y) Uh{y, rj, t)dydr]. 

We then deduce 



H^^ix,^, s) - HX\x,^,s) = {I- e-^^-)Tr{W, o + ... 

... + / W{x-y) (uh{y,v,t)-Vh{y,v,t))dydr]. 
In view of the preceding points, 

- W/l(.,t)||il(]R2n) < C / ||Vw?i(.,s)}||l,l(]R2n) /l + || U/i (. , s) - W/, (. , s) || 1,1 (]R2n ) 

From lemma 7.1 (with A; = 1) we then deduce (7.13). 

End of the proof of theorem 2.1. First step. With the above notations, (TDHF) equation is written as: 



ds. 



We consequently have 



ih^-^ = [Hr{t),P.{t)]. 



Since the propagator associated to this equation preserves the trace norm, and since 3^(0) = Th{ph{0)), we 
then deduce 

\\Ph{t) - B^mciin) < WPhiO) - ThiPhmWcH-H)- 
Consequently, from the proposition 4.2 (with m = 0) and theorem 5.1 (point c)), 

(7.14) {2nh)-^ar''[ph{t) - S^t)) ||li(r2„) < ||p^0) - ThiPhmWcHn) < CVhltiphiO)). 

Second step. We shall bound in norm Bfi{t) — Chit) and in this aim we shall show that Bh{t) and Ch{t) 
verify similar equations. The operator Bh{t) verifies (7.6) whereas 

ih'^^ = ih{2^hrOPr\{H^^{.,t),WH{.,t)})- 
With the notations of section 3 in [AKN] we have: 



Consequently, 



(7.15) i/i^^ - [Hr{t),CH{t)] = {2nhrR'j^\Hr{.,t),WH{.,t)). 

We know that R^^\f,G) = for all function G when F(x,^) = l^p. We can then replace Hf^^{.,t) by 
Vq{.,ph.{t)) in the right hand-side of (7.15). By combining (7.15) with equation (7.6) verified by Bh{t) and 
using Duhamel's principle and (7.4), we obtain 

1- I' "^5(2)/ 



\\Bh{t)-Ch{t)\\c^(n) < ^(W / \\Ri:>{V,{.,ph{s)),Wh{.,smc^H)ds. 
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From theorem 3.1 of [AKN] applied with N = 2, F = Vq{., ph{s)), G = Wh{-, s), p = oo, q = 1, we have 



a+^<6n+S 

:<>2,/3>2 



Since the potentials V and W are in M^°°°°(]R") and since the L^(]R^") norm of Wh{-, s) is bounded (lemma 
7.1) then the derivatives of all order of Vq{., ph{s)) are bounded. For ,9 > 1 the function W^Wh{-, s) verifies 
the estimates (7.11) of lemma 7.1. Consequently, when h G (0, 1], 



l|i?r(n(.pU5)).^/.(.s))IUM«)<cW^'^'^rU(o)). 



Consequently, 



\\Bhit) - C^micHn) < mVhlf^ipniO)). 
We then deduce (from the Proposition 4.2 with m = 0) that: 

< c{t)Vhii'{phm. 



Z,i(IR2'») 



(7.16) (27r/i)-" ar^'^iBhit) - Ch{t)) 

Third step. From the proposition 4.4, we have 

{2nh)-^arHCh{t) - £'ft(t))|Ui(iR-) < We^^'ivhi^t) - Wh{.,t)\\L^TR^-) < \\{vk{.,t) - wn{.,t)\\LHTR^-y 
Then, from lemma 7.2, 



(7.17) (27r/i)-"||ar'^(C^(t) - 2?/.(t))|Ui(R2„) < C7^(p^(0)) 



h+ ||wft(.,s) - w^(.,s)||i,i(]R2„)ds 



Fourth step. From proposition 4.4 we have 

(27r/i)-Vr*=(D^(i)) = {27rh)-^e^^aZ^'{DH{t)) = ei^Vh{.,t). 
Consequently, from the lemma 7.1 with k = 2, we obtain 

||^;^(.,i) - (27r/i)-"ar'=(i?ft(t)|Ui(H2„) < IKe*^ - /)t;ft(.,t)|Ui(R2„) < ... 



(7.18) 



... < C/i||VV(.,t)IUi(iR-) < Cit)y/hlZ{phm. 



Since Uh{-,t) is defined by (2.3), we obtain from estimates (7.14), (7.16), (7.17) and (7.18) obtained in the 
four steps that: 



\\uh{.,t)-Vh{.,t)\\LHM^-) < C{t)lf[{phm 



/ ||wft(-,s) - v/»(.,s)||ii(]R2,.)c;s 

JlO,t] 



From Gronwall's lemma, we have 

\\uH{.,t) - Vh{.,t)\\L^^M^^) < C(t)y)^/r(p^(0))e^(*)^''('"-(°)) 
a h G (0, 1], with a different constant C{t). Theorem 2.1 is complete. 



□ 



8. Approximative composition of symbols with the Wick calculus. 
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Being given two continuous operators A and B in iS(lR") we can define the Wick symbols crj^'*'^'^ ( A) , a]"^'^'' (B) 
and a'^^'^'^{A o B). When one of these two symbols a'j^'"^''{A) or crj^*'^'^ (S) is a polynomial function, we can 
express a'^^'^'^{A o B) using an exact formula with the following notations. For any function F e C^(IR^") 
we set 



OF .dF 



djF 



1 



dF .dF 
dxj d^j 



For any multi-index (a, 13) we set 

When F is a function of several variables in IR^" denoted by W, Y, etc., we shall write as a subscript the 
letter giving the variable on which act the operator. 

With these notations, if one of the two symbols a^^''''{A) or cr^*'=*(B) is a polynomial function then we have: 

(2/i)" 



U) 



ar>'{AoB){X)^ ^ 



d'^ar'iAhXX) 



d"ar\Bu){X) 



|a|<m 



This formula justifies the terminology since it carries the Wick order of creation and annihilation operators. 

In this section, we shall give a similar formula when one of the two operators is trace class while the other 

one associated through the Weyl calculus has a symbol in VK°°°°(IR^"). In this case, for all integer numbers 
m > 1 and for all ft, > we shall denote by Rm{A, B;.;h) the function on IR^" defined by the equality 

(2/i)" 



|a|<m 



(8.2) af^'=(^oB)(X)= V ^ 9Vf^'=(A)(X) d" ar''\Bh){X) + R^{A,B-X-h). 



The main result of this section is the following one. 

Theorem 8.1. We consider two operators A and B such that 

1. The operator A is written as A = Op1''y\F) where F e VF°°°°(IR2»). 

2. The operator B is trace class. 

Then, for all integer numbers m > 1 there exists a constant Cm > ( depending only on m and n ) such that 
the functions Rm{A, B; X; h) defined by the equality (8.2) and the analogous function Rm{B, A; X; h) verify: 



(8.: 



.3) (27r/i)-" / Rm{A,B-X-h) dX<C^n \\Bh\W(H) /i'"'/'||a"i^lUoo(K. 



where Emn is the following set of multi-indexes, 

(8.4) Emu = {a e IN", m<\a\< sup(m, n + 1)}, 
and similarly for Rm{B, A; X; h). 

We shall use the function Sh{A) defined in (4.2) and those which are similarly associated to B and Ao B. 
We shall refer as Sh{A) for the bi-Wick symbol of A. Using (4.5) we have: 

(8.5) ar^'=(Aoi?)(X) = (27r/i)-" / Bh{X,Y,Y,X) Sh{B){X,Y) Sh{A){Y,X) dY. 

The proof of theorem 8.1 relies on a Taylor expansion for the bi-Wick symbol of the operator Ah = Op^^^''{F) 
in a neighborhood of the diagonal. It is standard that this bi-symbol is holomorphic in X = a; + and 
anti-holomorphic in Y. We can give a quick proof. Iterating (4.5) we see that 



(8.6) 



Sh{A){X, Y) = (27r/i)-" / Bh{X, U, V, Y)Sh{A){U, V)dUdV, 
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where Bh is defined in (4.8). Replacing in (4.8) the expression (4.3) giving the scalar product of two coherent 
states, we obtain an explicit expression of the function Bh{X, U, V, Y), showing that it is holomorphic in X, 
anti-holomorphic in Y and Sh{A) inherits therefore of these properties. Prom these points, 



The first step in the proof of theorem 8.1 is an uniform estimation as h tends to of the remaining terms of 
order m in the above expansions. Namely, 

(8.7) Sh{A){X,Y)= ^ i^^^)l(d-ar'''^{A)){Y) + Rrn{A;X,Y;h), 

\a\<m 

(8.8) Sh{A){X,Y)= J2 ~sP^ (B^arHmX) + Rm{A; X, Y; h). 

We shall restrict ourselves to the case where A = Op^'^^^F) is a pseudo-differential operator. 

Proposition 8.2. Let A be an operator of the form A = Op^^"^'' (F) where F G W°°°° {Wl^"^) . For all integer 
number m > 1, let Rm{A;X,Y;h) and Rmh{A; X,Y;h) be the functions defined in (8.7) and (8.8). Then, 
there exists a function Gm in //^(M^") such that, for all h in (0, 1] 



(8.9) 



where Emn is defined in (8.4). An analogous estimate remains valid for the function Rm{A;X,Y;h) defined 
in (8.8). 

Proof. Proposition 4.4 gives an expression of the Wick symbol of Ah = Op^^'"\F) 

ar'^iAhKX) = (tt/i)-" / e-^^^F{Z)dZ 

Since Sh{Ah){X,Y) is a function holomorphic in X, anti-holomorphic in Y and equaling to cr'^^'^'' (Ah) on 
the diagonal, we necessarily have 

Sh{Ah){X,Y) = {Trh)-^ [ e-*(^-^) (^-^)j^(Z)dZ. 

In order to derive the Taylor expansion, we set for all 6 G [0, 1] 

Fh{e, X, Y) = (tt/i)-" / F{Z)ei'^^^'^'^'^Uz, 

ip{e, X,Y,Z) = -(Z -Y - 9{X - y)) .(Z - Y)dZ. 

Thus, we have Fh{l,X,Y) = Sh{Ah){X,Y) and Fh{0,X,Y) = Sh{Ah){Y,Y). The Taylor expansion is then 
written as 

Sh{Ah){X,Y)= J2 ^d^Fh{0,X,Y) + Rm{Ah;X,Y;h), 
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'■I r-t n\m—l 

Rm{Ah; X, Y-h)= / \ d'o"Fh{0, X, Y) dO. 



Differentiating the above equality through the integral and integrating by parts, we obtain: 
d^Fhie, X, Y) = (tt/i)-" / ei^^^'^'^'^^ \{X - Y).dz\ F{Z) dZ. 



Consequently, 



^dlF,{e,X,Y) = V P\ nh)- I e*^(^'^'^'^) (d'^F){Z) dZ. 



|a|=fe 



In particular, for = 0, 

(tt/i)-" / e-^*'(°'^'^'^' rfZ= (tt/i)-" / e^l^"^!' (d'^F\{Z) dZ 

The Taylor remaining term is given by: 

Rm{Ah; X,Y;h) = m V ^^~P" {nhy [ (1 - 0)m-i^j^^{e,x,Y,z) gap^^) dOdZ. 

i[0,l]xIR- 

Consequently, 



\Rm{Ah;X,Y;h)\ < m\X -Y\ 



We see that: 



\a\=m 



(nh)-- [ 
J[o, 



(1 _ 0^m-l^iRe^(e,X,Y,Z) ^q^^. 



Re^{e, X, Y,Z) = - 
If m > n + 1 and C > the following function 



\X-Y\ 



Gm{X) = C \X\ 



Jo 



e * 



is in L^(IR^"). We can find (7 > such that (8.9) is satisfied. When m < n we have 

i<|a|<n 

When m < n we can then find C > such that the inequaUty (8.9) is verified when setting 



,{Ah;X,Y;h) = Rn+,{Ah;X,Y;h)+ ^ a!^^" (aic7r"(A)) (F). 



G„(X) = G„+i(X) + Ce-^l^l' El^l 



This function Gm is also in L^(IR^"). In both of the two cases we have the estimation (8.9). The estimation 
concerning Rm{Ah; X,Y;h) is similarly proved. 

□ 
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End of the proof of theorem 8.1. Using (8.5) and using the asymptotic expansion in the first variable (here 
Y) of the function Sh{Ah){Y,X) given in the proposition 8.2 when = Op^^^\F), we obtain: 



ar'\AoB){X) = ^^ 9"ar°fc(A)(X) + R^iA, B;X;h) 

\a\<m 



where 



a^{X, h) = (27r/i)-" / Bh{X, Y, Y, X) Sh{B){X, Y){Y - X^dY 

R^{A, B; X; h) = (27rft)-" / Bh{X, Y, Y, X) Sh{B){X, Y)R^{A; X, Y; h) dY. 
Prom (4.9) we have: 

Bh{X,Y,Y,X) = Bh{X,X,Y,X) = e-^(^-^)(^-^). 

Therefore we obtain that 

a„(X,/i) = (2/i)l«l(27r/i)-" / Sh{B){X,Y)dxBh{X,Y,Y, X) dY 

= (2/i)l«l(27r/i)-"9j / Bh{X,X,Y,X)Sh{B){X,Y) dY 

= (2/i)l"l 9" (S)(X). 

We have used the fact that dxSh{B){X,Y) = and the fact that the bi-symbol Sh{B) verifies the property 
(8.6). Besides, from the bound of Rm{A; X, Y; h) given by proposition 8.2, 

\Rm{A,B-X-h)\<Cm{F,h) (27r/i)-"^^^ g-^l^-^l' |5^B)(X, F)|G„ (^^^^) dF, 

Cm{F,h)= /i'"'/'l|5"i^lUoo, 
where Gm is a function in L^(]R^"). Consequently, 



(27r/i)-" / Rra{A,B;X-h) dX <... 

... < CUF,h) (27r/i)-2"^^ e-Al^-^l' \Sh{B){X,Y)\Gm (^^) ^^^Y. 

From proposition 4.3, the point concerning AoB \s proved and the one concerning So A is similarly derived. 

□ 

9. The equation satisfied by the Wick symbol of a solution to (TDHF). 

If one only supposes that {t) is a classical solution to (TDHF) without assuming any additional hypothesis, 
then the Weyl symbol of ph{t) is a continuous function on M^" with derivatives understood in the sense of 
distributions. However, from proposition 4.2 it is noted that its Wick symbol is in W^°°^(]R^") and therefore 
it is C°° . We shall prove in this section that the Wick symbol satisfies a differential equation, with an 
asymptotic expansion in powers of h, the first term of the equation being the Vlasov equation, and the error 
term being as small as wanted, in terms of power of h, and of i^(IR^"). For the introduction ofsuch an 
equation, see [DLERS]. 



22 



Proposition 9.1. a) Let ph{t) he a classical solution to (TDHF). For alltGJR set: 
(9.1) ^,{X,t) = {27rh)- ar'\pk{t)) 



(9.2) 



Then, 



»„(x,j) = (2„ft)-"ar"(m(*)) 



a) We have, in sense of distributions on H^" x IR, 
(9.3) 

b) We have: 



We observe that when the potentials are turned off (V^ = = 0) then the Weyl symbol of Ph{t) is the only 
one that exactly verifies Vlasov equation (but in the sense of distributions). 

Proof. Point a) Let ph{t) be a classical solution of (TDHF). Let Uh{-,t) be the function defined in (9.2). 
Let V be a C"^ function with compact support in M^" x IR. For all teTR, set Ah{t) = Op^'''"\ip{.,t)). The 
mapping is from IR into £^ {H) and it is continuous from IR into V. We shall represent the integral 



7ft = (27r/i)-" / Mx,^,t) 



dxd^dt. 



We know that: 



ar'{h[A,A,itmx,0=2zj2^,^i^,^,t). 

Consequently, we have from (3.2), 

Ih = j^Tr{ph{t) o K(i) - ih[/^,Ah{t)]))dt. 

Since the operators ph{t) and Ah(t) are in I?, wc have: 

Tr{pn{t)[A,An{t)]) = -Tr{[A,phit)]Af,{t)). 

Consequently, 



Using (3.2) again. 



-h = j^Tr({p'^{t) - ih[A,ph{t)]) o Ah{t)yt 
= ^J^Tr{[VMt)),pHit)]oAh{t))dt. 

-h = ^(27r/i)-" / a];:'yH[VMt)),Ph{t)])^{x,^,t)dxd^dt. 
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Point b) Let Ph{t) be a classical solution to (TDHF). Let Uh{-,t) and (ffi{-,t) be the functions defined in 
(9.1) and (9.2). We know (see proposition 4.4), that Uh{-,t) = ei'^(ph{;t)- Applying the operator e^'^ to 
the two hand-sides of (9.3) we obtain (9.4). □ 

Combining theorems 9.1 together with 8.1 and lemma 7.2 together with proposition 4.4 wc directly obtain 
an equation satisfied by the Wick symbol of a solution to (TDHF) without any supplementary hypothesis 
and with an arbitrary high order of accuracy. However, without supplementary hypothesis, this equation 
does not allow to obtain an asymptotic expansion of the Wick symbol. 

Theorem 9.2. Let V and W be two potentials in W°°°°{TR"). Let ph{t) be a classical solution to (TDHF) 
such that Ph{0) > 0. Let Uh{-,t) be the function defined in (2.3). For all functions f in L^{M?^), set: 



^hU){x) = (e4^F)(x) + / W{x- y)f{y, r,)dydr). 



Then the function Uh verifies, for all integer numbers m > 2; 



dt 



- E 



{2hY 



ih '—f ai 

l<\a\<m 



d'^MM.t)) 9 M;t) - d''wh{.,t) d 'i>h{uh{.,t))\ +R^{.,t,h) 
where the function Rm{.,t, h) is in Z/^(]R^") for all f e H. For all T > 0, there exists Cm{T) > such that: 

(27r/i)-" / Rm{X;h) dX <Cm{T)h'^-'^. 



10. The analogue of Ehrenfest's time. 

The below result shows that the limit in corollary 2.2 cannot be uniform on IR, even without potentials V 
and W entirely vanishing, since exchanging there the two limits is false. In other words, according to the 
terminology of [BR], [dBR], the analogue of the Ehrenfest time for coUary 2.2 is finite. 

Theorem 10.1. Let h> be fixed. Let Ph{t) be a classical solution to the (TDHF) equation corresponding 
to V = W = 0. We suppose that ph{0) > and assume that the trace of ph{Q) equals 1. Let Uh{.,t) and 
Vh{.,t) be the functions defined in section 2. Then we have: 

(10.1) lim ||«ft(.,i)-v^(.,i)|Ui(iR2n) =2. 

We recall here that the two functions Uh{.,t) and Vh{.,t) are in the unit ball of X^(]R^"). 
Proof. When V = W = the function Uh{., 0) verifies from the proposition 9.1, 



(10.2) 




The function Vh{.,t) is the solution to the Vlasov equation equaling to u/j(.,0) for t = 0. Then, Vh{x,^,t) = 
Uh{x — 2t^, ^, 0) in the case of vanishing potentials. Set 

(10.3) Uh{x,t)= [ Uh{x + 2t^,^,t)d^. 
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We have: 

\\uh{;t)-Vh{.,t)\\LiiM^^)= I \uh{x + 2t^,^,t)-Uh{x + 2t^,^,t)\dxd^> I \Uh{x,t) - Uh{x,0)\dx. 
Prom (10.2) and (10.3), 

-g^{x, t) = -h j^^ J2 a^^^ + ^' = "^htA^Uhix, t). 

Consequently, 

UH{.,t) = e''*"^-Uh{.,0). 

It is standard that: 

hm We^^-F - F||ii(]R„) = 2||^^|Ui(m.). 

A— >+oo 

if a function f > is in L^(]R"). The proof of the theorem is therefore completed. 

□ 
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